Comment on "Absence of spin liquid in non-frustrated correlated systems" 
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In a recent Letter, Hassan and Senechal [1] discussed 
the existence of a spin-liquid phase of the half-fihed Hub- 
bard model on the honeycomb lattice. Using schemes, 
such as the variational cluster approximation (VGA) and 
the cluster dynamical mean field theory (CDMFT) in 
combination with exact diagonalization (ED), they ar- 
gued that a single bath orbital per site of the six-atom 
unit cell is insufficient and leads to the erroneous conclu- 
sion that the system is gapped for all nonzero values of 
the onsite Coulomb interaction U. In contrast, we point 
out here that, in the case of the honeycomb lattice, six 
bath levels per six-site unit cell are perfectly adequate for 
the description of short-range correlations. Instead, we 
demonstrate that it is the violation of long-range trans- 
lation symmetry inherent in CDMFT-like schemes which 
opens a gap at Dirac points. The gap found at small 
U therefore does not correspond to a Mott gap. As a 
result, present CDMFT schemes are not suitable for the 
identification of a spin-liquid phase on the honeycomb 
lattice. 

As shown in Ref. [2], the cluster self-energy ob- 
tained within ED CDMFT [3] using six bath levels is 
in nearly quantitative agreement with results derived 
within continuous-time quantum Monte Carlo (QMC) 
CDMFT [4]. As the bath in QMC is infinite, the self- 
energy is not subject to finite-size effects. The reason 
for the good agreement is that, because of the semi- 
metallic nature of the honeycomb lattice, the projection 
of the infinite-lattice bath Green's function onto a finite- 
cluster Anderson Green's function is not plagued by the 
low-energy-low-temperature disparities that arise in the 
case of correlated metals and the square-lattice Hub- 
bard model. Moreover, because of the hexagonal sym- 
metry, this projection can be performed in the diagonal 
molecular-orbital basis with non-symmetrical density of 
states components [3]. The main issue in Ref. [1] con- 
cerning the symmetry of bath levels then does not arise 
and the two independent bath Green's function compo- 
nents are fitted accurately using a total of six parameters 
(see Fig. 24 of Ref. [2]), while in the site basis only one 
fit parameter is available. 

Within CDMFT, the self-energy at Dirac points K 
exhibits the low-energy behavior [3]: 'Ej{K,iuJn) ~ 
aiuJn + 6^/[?w„(l — a)], where a is the initial slope of 
Im [Sii(ici;„) — Ei3(iw„)] and b the low-energy limit of 
He [T^ 12 (iuJn) — Si4(iw„)]. Here, is the self-energy 
in the site basis and a;„ are Matsubara frequencies. 
The excitation gap at temperatures T is given 
by A « 2|&|/\/l — a (see also Ref. [5]). Using the site 
notation: ai = (0,0), aa = (0,1), ag = (^3/2, 3/2), 
a4 = (V3,l), as = (V3,0), and ag = (^3/2, -1/2), 
the elements E12 and E14 within this cell are indepen- 


dent, so that A 7^ 0. However, site 1 is also connected to 
site 4 at (— \/3/2, —1/2) in the neighboring cell, requiring 
S12 = Ei4 which is not obeyed in CDMFT. Therefore, 
this violation of translation symmetry is responsible for 
the insulating contribution ~ 1 / iujn to the self-energy at 
K. Clearly, this term is not caused by the finite size nor 
the symmetry properties of the bath employed in ED. In 
fact, in view of the good agreement with ED, the self- 
energy in QMC CDMFT [4] exhibits the same insulating 
contribution, so that the density of states should also 
reveal a gap at small U and low T. 

To restore translation symmetry, we have performed 
ED calculations within the dynamical cluster approxi- 
mation (DGA). This approach ensures S12 = E14, giving 
semi-metallic rather than insulating behavior at small U 
[6]. The condition E12 = E14, however, cannot be gen- 
erally correct for correlations within the unit cell. The 
critical Coulomb interaction t/c ~ 6 found in DGA should 
therefore be used with caution. 

For interactions in the range J7 w 4 . . . 5, ED and QMC 
CDMFT [3,4] yield excitation gaps which agree remark- 
ably well with the charge gap derived in large-scale QMC 
calculations [7]. Evidently, this gap corresponds to a 
Mott gap induced by short-range correlations and is only 
weakly affected by the lack of long-range translation sym- 
metry. On the other hand, the key question as to how this 
gap closes, i.e., how the semi-metallic phase with weakly 
distorted Dirac cones is recovered at small U, cannot be 
studied adequately within CDMFT. The gap obtained 
at small U is an artifact caused by the violation of long- 
range lattice symmetry and does not represent a true 
Mott gap. As a result, the identification of a spin-liquid 
phase on the honeycomb lattice is presently not feasible 
within CDMFT-like methods. 

In conclusion, the origin of the narrow excitation gap at 
small U found in ED CDMFT for the honeycomb lattice 
is not the finite bath as stated in Ref. [1] but the lack 
of translation symmetry. This problem is also shared by 
QMC CDMFT and can be avoided, for instance, within 
ED/QMG DGA. 
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